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ft. A restricted dynamics, previously introduced in a kinetic model for relaxation phenomena in 

linear polymer chains, is used to study the dynamic critical exponent of one-dimensional Ising 
models. Both the alternating isotopic chain and the alternating-bond chain are considered. In 
^^ ■ contrast with what occurs for the Glauber dynamics, in these two models the dynamic critical 

exponent turns out to be the same. The alternating isotopic chain with the restricted dynamics 
is shown to lead to Nagel scaling for temperatures above some critical value. Further support is 
given relating the Nagel scaling to the existence of multiple (simultaneous) relaxation processes, the 
dynamics apparently not playing the most important role in determining such scaling. 
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I. INTRODUCTION. 



The scaling hypothesis of Halperin and Hohenberg jD relates the time scale r and the correlation length £ and 
introduces the dynamic critical exponent z. In the case of Ising models, z was for a long time believed to be universal, 
depending on the nature of conserved quantities and of those features, for instance dimensionality 0, which determine 
their static universality class. However, it is now well established that for some simple systems this exponent is non 
universal H]- M. In particular, in the case of one-dimensional Glauber dynamics EQ| the alternating isotopic chain 
q , J7J presents universal behavior (in the sense that it leads to the same value of the dynamic critical exponent as the 
homogeneous chain) whereas the alternating-bond chain does not m- p8[ (see however Ref. [}L9|), This breakdown 
of the dynamic scaling hypothesis is due to the different energy barriers determining the diffusion constant and the 
correlation length and hence entering into the value of the dynamic exponent z. This exponent is also typically 



strongly affected when the dynamics takes place on a self-similar background. For instance, the study of the Glauber 
dynamics in branching and non-branching Koch curves as well as in the Sierpinsky gasket [ p0| indicates that z is non 
universal. Nevertheless, an analysis of the critical Glauber dynamics on the Fibonacci-chain quasicrystal || shows 
that the dynamic exponent z is identical to that obtained for the alternating-bond Glauber chain. Another interesting 
line of research that has also been followed M is to understand how different dynamics affect the universality. This 
work is partly aimed as another contribution to such an understanding. 

Kinetic Ising models have also been used in a variety of contexts not limited to critical dynamics. In fact, in a 
previous paper |2lf| we introduced a quasi one-dimensional kinetic Ising-like model to study relaxation phenomena in 
linear polymeric chains, including a region close to the glass transition. In our original model the chains are made up 
of N segments, each of which may be found in two possible orientations and the Hamiltonian was chosen as to reduce 
to the one giving the intramolecular energy of the Gibbs-di Marzio lattice model p2f . For the stochastic dynamics, 
we adopted a rule of transition for the configurational changes which was tied to the creation or disappearance of 
flexes. As a consequence of this restricted dynamics, only some states are selected, and in the magnetic language, this 
implies that in the model the domain wall motion is through a biased random walk. The purpose of this paper is, 
on the one hand, to introduce the restricted dynamics into the two one-dimensional alternating kinetic Ising models 
that have exact solutions with the Glauber dynamics, namely the isotopic alternating chain and the alternating-bond 
chain. On the other hand, we will use these extensions to examine the dynamical critical exponent and, in the case 
of the isotopic chain, whether the inclusion of several relaxation times may be related to multifractal behavior. This 



last goal is motivated by the fact that recent dielectric susceptibility measurements in several glass-forming systems 



and covering wide ranges of temperature and frequency |23]- 25 1, have suggested that the master curve in which all 
measurements are shown to scale (Nagel plot) may be understood in terms of multiple relaxation processes. In fact 
the existence of simultaneous relaxation processes is assumed to be connected to multifractality much in the same 
way that this concept is present in theories of chaos. Very recently we have shown that the alternating isotopic chain 
with Glauber dynamics leads to Nagel scaling P6| and the question arises as to the role that a particular dynamics 
may play in such a scaling. By examining the model with the restricted dynamics, we aim at shedding some light 
into this issue. 

The paper is organized as follows. Section 2 deals with the alternating isotopic chain while section 3 is devoted to 
the alternating-bond chain. We derive the wave-vector frequency dependent susceptibility the equations governing 
the time evolution of the wave-vector dependent magnetization and the critical exponent z in both cases. In section 
4 we address the issue of whether the restricted dynamics leads to the Debye-like scaling of the response functions (as 
in the uniform chain) or whether the so-called Nagel plot is useful in the case of the isotopic alternating chain. We 
close the paper in section 5 with some concluding remarks. 

II. THE ALTERNATING ISOTOPIC CHAIN. 

The model consists of a closed linear chain with N sites occupied by two isotopes (characterized by two different 
spin relaxation times) that are alternately arranged. The Hamiltonian is the usual Ising Hamiltonian given by 

N 

H= -J^Vj-o-j+i (1) 

where tjj is a stochastic (time-dependent) spin variable assuming the values ±1 and J the coupling constant. The 
configuration of the chain is specified by the set of values {<xi,c"2, ...crjv} at time t. As it will be argued below, the 
Hamiltonian will turn out to be not all that relevant in the specific calculations. Nevertheless it is necessary to define 
the states involved in allowed transitions. 

Instead of considering the Glauber dynamics, we assume a kind of transition associated to the motion of domain 
walls. The idea is similar but not identical to previous work by others in which either the domain wall motion is 
strongly suppressed at low temperatures || [Q or it is through a one-dimensional random walk (27]]. In our case the 
transition associated with the ith spin takes the form 

Ti {fTl, ...,<Ti,<T i+ l, ...,(Tjv-i,o'jv} — ► 

{d, ...,cjj_i, —<Xi,<Xi+i, ...,ctn-i,&n} , (2) 

and we impose a biased random walk for the domain wall motion |2q ]. In order to also account for the presence of 
the isotopes, we take the transition probabilities to be given by 

Wi(<7i-i,tJi) = (Xi(l - 7<7t_l<7i). (3) 

Here, 7 = tanh I -^-^ ) , ks being the Boltzmann constant and T the absolute temperature, and on is the inverse of 

the relaxation time Tj of spin i in the absence of spin interactions. It should be pointed out that the rule of transition 
stated in Eq. (H)allows for single site excitations such that the transitions are not correlated. Further, the choice 
made in Eq. (||) immediately implies that detailed balance does not hold for this model and also that not every state 
in the phase space of the system is accessible. In particular, one should note that within this model no equilibrium 
state exists, although when the time goes to infinity a steady state is eventually attained. This feature is shared by 
other purely stochastic models such as the biased Ising lattice gas and its many variations 129]. 

If we now let a\ and «2 represent the inverses of the free spin relaxation times of chains composed solely of spins 
of species 1 or species 2, respectively, then we can set on — 75 1 — (— l) l oi , where 5j = ai +" 2 and a 2 = ai ~ a2 . 

The time dependent probability P(<J\, 02, ■■■<Jn] t) = P ({&} , t) for a given spin configuration satisfies the master 
equation 

N 



dP({a N },t) ^, . , . . 
dt =-^2w i (v^ 1 ,* i )P({<T N },t) 



1=1 

N 



Y,w l (a l - 1 ,-<j l )P(T l {<j N } ,t) . (4) 



The dynamical properties we are interested in require the knowledge of some moments of the probability 
P ({c^} , tj. Hence, we introduce the following expectation values and correlation functions defined as: 

q i (t) = {a i (t))= 5>ii> ({*"}, t) (5) 

{**} 

r id (t) = {a t (t) <,&)) = £ a t ajP ({a N } , t) (6) 

{ CT « } 

and 

dj (*', t' + 1) = (t) (a t (t') g 3 (t' + t)) = 

53 ,:p(K'},t> 3P (K}iK'},t) (?) 

{<x«},{ CT «, } 

where O (t)is the Heaviside step function and the sums run over all possible configurations compatible with our 
rule of motion. The second equality of Eq. (pi), which gives the formal definition of the time-delayed correlation 
function, involves p ({tr"} | {c^'} ,i), the conditional probability of the chain having the configuration {c^} at time 
t' + 1 provided it had the configuration {<J N/ } = {a[, <j' 2 , ..., a-' N } at time t 1 . Multiplying the master equation by the 
appropriate quantities and performing the required summations we obtain the set of time evolution equations that 
will be used in our later development. These are given by 

^jf = -2«i (<& - 7fc-i) (8) 

and 

daj (t\t> + t) = ^ {f , )m _ 2 ^ ^ ( ^ t , + t) _ ^ ^ ( ^ fl + t) ^ (9) 

We now impose translational invariance and introduce §&, the (spatial) Fourier transform of qj, the t' — > oo limit 
of the (temporal) Fourier transform of q {if , t' + t) = a t j (£', t' + t) (with I = j — i) denoted by q(w), and Cfe(w), the 
spatial Fourier transform of cj(u;), defined through 

Qj = -7= 5Z ?fe exp (ikj) , (10) 

1 f°° 
c t (uj) = lim — / ci (f , t + 1) exp(— itot)dt (11) 

t'^oo 2tt J_ 00 

and 

C fc (w) = (cr_ fc cr fc ) w = — y^ ci(w) exp(-ikl). (12) 

l 

In terms of these quantities, Eqs. (0) and <M) may be rewritten, respectively, as 

^=MA (13) 

and 

ioJci(u) = rf - 2a t (cj(w) - yci-i(w)) , (14) 

where 



M* 



-2a 

2a 2 



i(l- 7 e- jfc ) 2a 2 (1 + je~ lk ) \ 
(1 - ie - lk ) -2ai (1 + 7e~^) / 



(16) 



limt_ >00 r;(t) is the value of the pair correlation function corresponding to the stationary solution of the equations 



of motion in the limit t — > oo. 



The solution to Eq. (13), which yields the magnetization, is straightforward, namely 



* fc (t) = e m ^ k (0) . 



(17) 



The relaxation process of the wave- vector dependent magnetization is determined by the eigenvalues of Mt . These 
are given by 



\i 



-2ai ± 2 



V^M^ 



(i 



f+Ll £ S3 £l%t\t \ 



(18) 



It should be stressed that the eigenvalues A^ contain both real and imaginary components. Therefore, the relaxation 
related to the real part of the eigenvalues will in general be modulated by the imaginary component, and care should 
be taken in defining an adequate correlation length for the critical dynamics. In fact, one can also associate this 
correlation length to the modulation of the oscillations. Such correlation length diverges at the critical point and 
hence the modulation eventually dissapears. This kind of behavior has already been thoroughly discussed in the 
context of the one-dimensional isotropic ferromagnetic XY model in an inhomogeneous transverse field pOfl . In any 
case, the inverses of the (fc-dependent ) relaxation times r^ of the ±k th modes are obtained from the real part of A^f. 



In the critical region, that is when T — > and k — > 0, A,. 



-4ai while \1 



0. This means that the critical mode 



is the one corresponding to Aj . As for the relaxation time, in this limit one gets 



Re(-A+)=-^ 



4(5?- 



7T72 yr2\ 



Ol 



m 2 



(19) 



where we have identified the correlation length £ as £ ~ e k b T by comparing the former expression with the one of 
the dynamic scaling hypothesis — ~ £,~ z f{£,k) . Therefore we find z = 2 , which is precisely the same result as 
for the alternating isotopic Glauber chain |7j. Note, however, that this correlation length corresponds to the one 
of an Ising model with an effective exchange constant i . This can be easily seen by noting that the equations of 
motion for the two-spin correlations are in our case formally identical to those in the Glauber chain, but in the latter 

7g = tanh I -^4p ) . Interestingly enough £ also corresponds, as expected, to the correlation length of the steady state 

attained by the system. The role of the effective constant ^ had been already pointed out in the case of the model 
introduced in ref. pifl . 

Now we turn to the calculation of the other interesting response function, namely the frequency and wave-vector 
dependent susceptibility Sk (w), which, by virtue of the fluctuation-dissipation theorem |31j|, is defined by 



Sk (w) 



\OkO-k) 

k B T 



iu; (cr fc cr_ fc ) 
knT 



(20) 



where (crfc<j_fc) =1/(1 — 7COsfc) cosh^py is the static correlation function and {akU-kj^ the Fourier transform of 
the dynamic one. After some rather lengthy but not too complicated algebraic manipulations starting with Eq.(nJ) 
we arrive at 



S k {oj) 



1 



ksT (1 — 7 cos k) cosh -^-^ 



- ik 



iui (iuj + 2ai(l 



7e 



)) 



{iu + 2ai) - Aal + 47 2 (a\ - a\) e 



(21) 



The result embodied in Eq. (plh , constitutes the proper framework in which to discuss the issue of Nagel scaling 
in the relaxation of a linear chain with translational invariance and restricted dynamics. This will be postponed until 
section 4. For the time being, we just quote the equivalent result for the isotopic chain with Glauber dynamics [R, 
p6fl, namely 



s fc G M 



k B T (1 - 7g cos k) cosh -^ 

iui {iuj + «i(l + 7g cosfc)) 
(iw + ai) - 57 G aia 2 (1 + cos2fc) — a| 



(22) 



which is also required for such a discussion. In the next section we will consider the other simple alternating Ising 
chain including the dynamics allowed by our rule of transition. 



III. THE ALTERNATING-BOND CHAIN. 



The model consists of a closed linear chain with N sites and it is characterized by different coupling constants Jj 
(to be specified below). N is taken to be an even integer and periodic conditions are also imposed. This model can 
be applied in the description of dimerized structures and the Hamiltonian is again of the Ising type, i.e. 



H = - 



N 

E 



JjVjVj+l 



(23) 



where the parameters and variables are defined as in the previous section. Once more and for the same reasons as 
above, we specify the Hamiltonian to have a precise definition of the states that intervene in allowed transitions, but 
it will play no further relevant role in the calculations that follow. 

Considering again the same kind of biased random walk for the domain wall motion and the rule of transition given 
by Eq. (0), in this case the transition probabilities are taken to be 



Wi((Ti-i,<Ti) = a(l - Pi-KTi-KTi) 



(24) 



where a is the inverse of the free spin relaxation time, 0j — tanh f -7-i— J and Jj = A ( J\ + J^j 

J\ and J2 represent the coupling constants of two different uniform Ising chains, respectively. It is also convenient to 



( -^- (Ji - J 2 ). Here 



introduce the quantities 71 = 5 tanh I -r^f ) + tanh 



J 2 

k B T 



and 72 = \ 



tanh 



k B T 



-tanhf-^ 



so that f3j may 



be expressed as 0j 
read 



and 



7i 



(— 1) J 72. The analogous forms of Eqs. (pi) and ( 0), obtained using a similar procedure, 



dqj . . 

— = -2a{q j -p j - 1 q j -i) 



dc ZJ (t\t' + t) = r ^ {t r )m _ 2a {cii {t/j t , +t) _ p._ lC . t ._ x (t ' ? t > + t)) . 



(25) 



(26) 



As in the previous case, we now impose translational invariance and again use g%, Cj(w), and C^{uS) defined as in 
section 2 to rewrite Eqs. (g5|) and ( Pq) as 



d&k 
dt 



N fc $ 



fcv/,. 



and 



where 



lUJCl(U>) 



2a(ci(oj) -j9j_iq_i(w)) 



(27) 



(28) 



$*: = 



qk 
qk-r 



(29) 



N fc 



-2a (1 



71 e 



— ik\ 



—ik 



2a^e- lk -2a (l ' ~ ■ ~-^ 



-2a72e 
71 e 



(30) 



The solution to Eq. (27) , which yields the magnetization, is again straightforward, namely 

$ fc (t) = e Nfct <I> fc (0) . 
The eigenvalues of N^ may be obtained very easily with the result 



^k 



-2a I 1 ± \l-fle- 2lk - 7|e" 2lfc 



(31) 



(32) 



and again contain real and imaginary components. Proceeding as in the derivation of Eq.(n9h, we find that the critical 
mode corresponds to A fc and that the relaxation time is given by 



Re -A 



i ~ 2 «r 2 



1 



m 2 



(33) 



Here, considering that 3\ > J2, the correlation length £ has been identified as £ 
(p3J) with the dynamic scaling hypothesis J- ~ £ /(£&)• Therefore, we also find z 



■'2 



e^B 1 by again comparing Eq. 

2 in this case, which coincides 
with the result derived above for the isotopic chain and, consequently, the restricted dynamics produces a dynamic 
critical exponent which is independent of the kind of interactions in these models. Note that once more the correlation 
length £ corresponds to the one of an Ising model with an effective constant 4p. 

The dynamic scaling exponent z may also be derived in this case by an alternative argument based on domain wall 
motion. Consider the limit T — > in which critical slowing-down occurs and domains are formed. According to Eq. 
(|24j), if spin j is in the interior of the domain, in this limit Wj(aj—i, <Jj) — > irrespective of the value of /3j_i;on the 
other hand, if spin i belongs to a domain wall, u>j((7j_i,(7i) — ► 2a also independently of /3j_i. Therefore in the low 
temperature limit Wj << w%, so that the dynamics of the system is equivalent to a biased random motion of the 
domain wall for the Hamiltonian system, in this instance the Ising model. Hence, within restricted dynamics, the 
relaxation time To for the decay of a domain of size £ (which is the correlation length associated to the steady state 
of the stochastic model) is such that To ~ £ and may be related to a biased random walk [Q . To this end, let Nq be 
the number of spins that must flip to obtain a domain of size £j (the correlation length of the Ising model). Since the 
domain has the size of the biased random walk with Nq steps and z = 1 for the biased random walk, £j ~ Nq. On 



N w7 



£r/(2a). But £1 = £ , so that z — 2 for the alternating-bond model for any choice of 



the other hand, To 
J 1 and J2. 

Finally, for the sake of completeness we will also compute the frequency and wave-vector dependent susceptibility 
Xk ( CJ ) °f this model. Taking the same steps as in the case of the isotopic chain, we find 



S k (w) = 



{<7 k <7- k ) x 

k B T 



-ik\ 



1 



iu> (iu) + 2a + 2ajie' 
(iuj + 2a) 2 + 4a 2 e~ 2ifc (7! - 7?) 
2iuja r y2&~ lk 



k B T 



(iu + 2a) + Aa 2 e~ 2lk (7 



7L 



with 



CTfc(7_fc) = 



(1 - uiu 2 ) [1 + U1U2 + [u\ + U2) cosfc] 
1 + u\u2 — 2uiU2 cos 2k 



(34) 



(35) 



[VkV-k+TTJoo - 



i (1 + U1U2) (v,2 — ui) sin k 
1 + ufu 2 , — 2u\U2 cos 2k 



(36) 



and Uj = tanh f 



■h 



2k B T 



j = 1,2. Again for the sake of comparison and in order to correct a misprint in the formula 



that appeared in ref. |7J] , we quote the equivalent results using the Glauber dynamics, namely 



#M 



\ak0-k) 
knT 



\OkV-k+n) 



ioj ( iuj + a + 2«7 cos k J 



(iu + af - 2a 2 (( 7 ' G ) 2 - (,5 G ) 2 ) cos2fc - 2a 2 (( 7 ' G ) 2 + (8 G ) 2 ) 



2a8 G uu sin k 



k B T 



(iu + a)' - 2a 2 (( 7 ' G ) 2 - (5 G ) 2 ) cos2fc - 2a 2 (( 7 ' G ) 2 + (<5 G ) 2 ) 



(37) 



with 

,G _ (1 - tgtg) [1 + uf ^ + (uf + tg) COS k] 

[a k a_ k)oo - i + ( u G u G )2 „ 2u G u G cos2fc . ^»; 

G _ ^(l + M f M G)( U g- M f)sinfc 

^ fc CT_ fc+w ; oo _ 1 + (u G u G )2 _ 2u G u G cos2fc > W 

7'° = itanh ( J fc + j 2 ) ) ^ G = -^tanh ( J ^j 2 ) and u^ = tanh (x^r )> J = I) 2 - This concludes our analysis of the 
alternating-bond chain with the restricted dynamics. 

IV. SUSCEPTIBILITY AND SCALING BEHAVIOR. 

Thus far, we have examined the behavior of the critical dynamic exponent. In this section we will continue our 
exploration by considering the scaling properties of the susceptibility in the isotopic alternating chain. This is most 
conveniently done through the quantity x( w ) — ksTSo(uj)/ (coctq) . It should be noted that if we look at the 
uniform chain, i.e. we set a\ = a 2 in Eq. ( pl| ) and take the limit k = in the resulting expression, \ obeys the usual 
normalized Debye scaling : if one scales the frequency with the inverse of the (single) relaxation time and the real and 
imaginary parts are then divided by their values at zero and one, respectively, one gets a universal curve. However, 
as mentioned in the Introduction, recently experimental work on dielectric relaxation has been reported in terms of a 
new scaling function [£3|- |E5) which is thought to be related to multifractal scaling. In this new scaling, the abscissa 
is (1 + W) log 10 (uj/ujp) /W" and the ordinate is log 10 (x"(u))u> p /u)Ax) /W '. Here, x" is the imaginary part of x(w), 
W is the full width at half maximum of x" , u) is the frequency and uj p the one corresponding to the peak in x" 1 and 
A-x = x(0) — Xoo is the static susceptibility. It is therefore interesting to see whether the isotopic alternating chain 
with the restricted dynamics leads to Nagel scaling in the same way that the model with Glauber dynamics does [£6| . 

Using Eqs. (ttq) and (plh with k = 0, x can be expressed in the form 



(1 -7)(ai + a 2 ) 

xM = — 



l-/(ai,a 2 ,7) l + /(ai,a 2 ,7) 



2 
Here, the (temperature dependent) function /(aii,a 2 ,7) is given by 



— Aq ilu — A 



(40) 



ti \ (ai-a 2 ) 2 -4aiQ 2 7 2 

/(ai,a 2 ,7) = ^=. (41) 

(ai + a 2 )v (ai - a 2 ) 2 + 4«ia 2 7 2 

It should be pointed out that for the case 7 = 0, we get 

X 7 =o = . ,\ + . * , (42) 

vjj + 2ai iui + I0L2 

so that the general structure of the result for the susceptibility of the alternating isotopic chain is preserved irrespective 
of the value of 7 (i.e. of the temperature), namely a linear combination of two Debye-like terms. 

With the aid of Eqs. ( [40|) and ( ffl| ) in Figs. 1 to 4 we present Nagel plots for the cases a.\ = a 2 = 1, <X\ = 1 
and a 2 = 2, a\ — 1 and a 2 = 100, and a% — 1 and a 2 = 1000, respectively, and different values of 1/T* = J/fc^T. 
We also include in these figures plots of x"( w ) vs - w/wp which are the natural variables of the Debye relaxation. 
While the first case (which as stated above corresponds to Debye behavior) does not show Nagel scaling, the situation 
somewhat improves in the second one (where clearly improvement means less dispersion in the curves), and when 
the two relaxation times are not only different but very far apart (third and fourth cases) the scaling is virtually 
perfect, provided the temperature lies above some certain critical value. For comparison, in Fig. 5 we show parallel 
results computed for the alternating isotopic chain with Glauber dynamics (c.f. Eq. (P2|)), for the case a\ = 1 and 
a 2 = 100. The similarity of the results of both models provides support to the idea that, irrespective of the specific 
dynamics, the coexistence of different relaxation mechanisms lies behind the Nagel scaling and that this only occurs 
if a threshhold temperature is surpassed. Nevertheless, the character of the different dynamics manifests itself in that 
the slopes of the decaying parts of the curves in the Nagel plot differ and that, for the same values of a,\ and a 2 , the 
scaling is more closely followed by the model with the restricted dynamics (c.f. Figs. 3 and 5). 

It should be noted that, as the insets of these figures indicate, both types of dynamics lead to two peaks in x"( w ) 
due to the presence of two different relaxation times. This feature has been also observed experimentally by Dixon et al 
p3[ , Lesley-Pelecky and Birge |24| and Wu et al [£5| in different materials and associated to the a and (3 relaxations. 
Very recently it has also been confirmed in experiments by Brand et al |32j who have also noticed the presence of a 
third relaxation process unexplained so far. 



V. CONCLUDING REMARKS. 

In this paper we have addressed relaxation processes and the question of the value of the dynamical critical exponent 
in kinetic Ising models on alternating linear chains. Two different issues have been examined in this context. In the 
first one, we have shown that with the restricted dynamics implied by the transition probabilities (c./. Eqs. (0) and 
(E4)), both in the alternating-bond chain and in the isotopic chain the dynamic critical exponent z turns out to be 
exactly two. This does not occur if the Glauber dynamics is employed. Hence, the value of the dynamic critical 
exponent and the dynamics implied by the rule of transition in kinetic models are deeply related. As for the second 
issue, the analysis of the Nagel plots in the case of the alternating isotopic chain indicates that the presence of at least 
two different relaxation mechanisms is required for the scaling of the susceptibility. This feature agrees with what one 
finds with the usual Glauber dynamics |2q| , as well as the appearance of plateau regions in the plot if the relaxation 
times are widely separated and of the existence of a critical temperature below which the scaling is not followed. One 
may reasonably wonder at this stage whether the alternating-bond chain with the restricted dynamics also obeys the 
Nagel scaling. Since for high temperatures this model is equivalent to a chain with a single relaxation time, it is not 
surprising that the scaling is not followed in this instance. This we have confirmed numerically for a variety of values 
for J i and J 2- 

Although our results for the isotopic chain suggest that the dynamics seems not to play a key role for the scaling 
to hold, one should bear in mind that the isotopic chain with restricted dynamics containing only two isotopes 
is somewhat peculiar, so that no definite conclusions on this issue can be reached at this stage. Concerning the 
differences, it is conceivable that the stochasticity may well be behind the fact that the slope of the decaying parts of 
the curves in the Nagel plots is larger for the restricted dynamics than for the Glauber dynamics. In this respect, it is 
important to point out that the experiments in which the Nagel plots have been more succesful concern glass-forming 
systems and that our original model was built so as to reflect the topological constraints that are assumed to be 
crucial in such systems. It is clear that the restricted dynamics related to such constraints is also not enough to get 
the scaling, as exemplified by the case of the uniform chain oc\ = U2 = 1. The true alternating isotopic chain with 
restricted dynamics examined here, on the other hand, includes both ingredients and provides a bona fide microscopic 
model in which the Nagel scaling is shown to arise. Moreover, the fact that we get a bigger slope in the restricted model 
is consistent with the experimental finding that such a slope is bigger for the orientationally disordered crystalline 
phase of cyclo-octanol (c.f. Lesley-Pelecky and Birge p3) than in the originally studied linear polyme rs [p5 [ ] . 

Notwithstanding the limitations of this model, our expectation is that both our earlier results (2M , fl26fl as well as 
the present ones provide some insight into the physical origin and validity of the hypothesis concerning the need of 
the simultaneous presence of multiple relaxation mechanisms as related to the proposal of the Nagel plots. A future 
challenge is to examine whether the Nagel scaling is also present in other Ising models recently studied in connection 
with glassy dynamics [B3[ which are based on the spin facilitated models originally introduced by Fredrickson and 
Andersen [[34|. Finally, one can conjecture that the appearance of a third relaxation process as observed in the recent 
experiments of dielectric relaxation by Brand et al p2 | may be hopefully catered for within our model through the 
inclusion of a third relaxation time. The investigation of this conjecture is presently in progress. 
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Figure captions 

Figure 1. Nagel plot for the case of a uniform chain («i — a^ — 1). Note that the proposed scaling does not hold 
in this case while, as the inset shows, the susceptibility obeys the usual Debye scaling. 

Figure 2. Nagel plot for a\ — 1 and cxi — 2. Here one can see an improvement of the scaling behavior as compared 
to the uniform chain case, except at low T*. In the inset the plot to test the performance with respect to the Debye 
scaling is presented. 

Figure 3. Same as Fig. 2 but for the choice ot\ = 1 and a2 = 100. Except at low T* values, the trend of improvement 
of the agreement with the Nagel scaling is apparent, while the opposite happens with respect to the Debye scaling. 

Figure 4. Same as Figs. 2 and 3 but for ot\ = 1 and a-2 = 1000. The scaling is virtually perfect in this case, except 
at low T* . Notice the explicit appearance of a plateau region in the plot. The behavior is here definitely non-Debye. 

Figure 5. Nagel plot for an alternating isotopic chain with Glauber dynamics with ot\ = 1 and 012 = 100 and different 
values of the parameter Tq = IzbT/2J. Here, x G (w) = kBTSQ(oj)/(aoao) 00 with (cqCo)^ = 1/(1 — 7G)cosh ^r 



and Sq(uj) computed with the aid of eq. (E3). Note the similarity of these results with respect to those of Fig. 3. 
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